Thouless numbers for few-particle systems 
with disorder and interactions 
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Considering N spinless Fermions in a random potential, we study how a short range pairwise 
interaction delocalizes the iV-body states in the basis of the one-particle Slater determinants, and 
the spectral rigidity of the iV-body spectrum. The maximum number grjv of consecutive levels 
exhibiting the universal Wigner-Dyson rigidity (the Thouless number) is given as a function of 
the strength U of the interaction for the bulk of the spectrum. In the dilute limit, one finds two 
thresholds: When U < U c i, there is a perturbative mixing between a few Slater determinants (Rabi 
oscillations) and g N oc \U\ P < 1, where P = N/2 (even N) or (jV + l)/2 (odd N). When U = U a i, 
the matrix element of a Slater determinant to the "first generation" of determinants directly coupled 
to it by the interaction is of the order of the level spacing of the latter determinants, gjv ~ 1 and 
the level spacing distribution exhibits a crossover from Poisson to Wigner, related to the crossover 
between weak perturbative mixing and effective golden-rule decay. Moreover, we show that the 
same U c \ signifies also the breakdown of the perturbation theory in U. For U c i < U < U C 2, the 
states are extended over the energetically nearby Slater determinants with a non-ergodic hierarchical 
structure related to the sparse form of the Hamiltonian. Above a second threshold U C 2 , the sparsity 
becomes irrelevant, and the states are extended more or less ergodically over <?jv consecutive Slater 
determinants. A self-consistent argument gives gN oc compare our predictions 

to a numerical study of three spinless Fermions in a disordered cubic lattice. Implications for the 
interaction-induced JV-particle derealization in real space are discussed. The applicability of Fermi's 
golden rule for decay in this dilute gas of " real" particles is compared with the one characterizing a 
finite-density Fermi gas. The latter is related to the recently suggested Anderson transition in Fock 
space. 



PACS numbers: 72.15, 73.20 
I. INTRODUCTION 

For non-interacting electrons the Thouless energy E c 
has proven to be a very relevant energy scale for sev- 
eral physical properties. The related "Thouless number" 
gi = E c / Ai, where Ai is the single-particle level spacing 
at the Fermi energy, plays an important role. For the 
disordered case, E c — fiD/L? , D being the diffusion con- 
stant of the electrons and L the relevant sample length. 
In this case, the Thouless number g\ is equal! to the di- 
mensionless conductance, i.e. the conductance in units 
of e 2 /%. This important relationship is at the basis of the 
scaling theory of localizations, which has been quite suc- 
cessful in describing transport in disordered metals for 
non-interacting electrons. 

E c is also an important energy scale for the spectral 
correlations of diffusive non-interacting particles in a rang 
dom potential. It was found by Altshuler and Shklovskiitl 
that the usual random-matrix correlationsu of the den- 
sity of states at different energies E and E' hold only 
when the relative energy \E — E'\ is smaller than E c . 
This means that one has only g± consecutive one-particle 
levels which exhibit the universal Wigner-Dyson rigidity. 
For \E — E'\ > E c a novel spectral correlation function 
was obtained which depends on the dimensionality and 
the diffusion constant. This new dependence rf-nd the 
crossover associated with it follow rather easilyH from a 



semiclassical theory for the spectral correlations due to 
Berryfj. _ 

When electron-electron interactions are introduced!!], a 
single-electron (or hole) excitation with an energy e ac- 
quires a finite width r sp (e)Q This width obviously in- 
creases with e. It has been calculated for au_ isolated 
metallic (gi ^> 1) compact quantum dot in Ref£3, where 
it turns out that at the Thouless energy E c this width 
becomes comparable to Ai and the single-quasi-particle 
excitations can no longer be resolved. Thus the number 
of single-particle levels that can be resolved is of the or- 
der of the Thouless number g\ . This result which agrees 
with the experimental findings of RefHj is universal and 



*The width of the single-quasi-particle excitation exists only 
when the energy of the excitation is high enough, namely 
above a threshold e* , and the golden-rule formulation for 
decay is valid. An analysis of this question either by using 
the formulation of Refill, p« via a correspondence with lo- 
calization on a Cay ley treeU, shows that this crossover en- 
ergy e* is of the order of %/E c Ai ~ Ai^/fTf. Note that since 
%/ E c Ai <§C E c , the level width is well defined and valid on the 
scale E c and for a large window below it. In Ref.cil further 
extremely interesting results were obtained to which we shall 
return later. 
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does not depend on material parameters, nor on the di- 
mensionality of the dot. 

The problem of interacting particles in a random po- 
tential is otgreat fundamental interest. In particular, the 
suggestion^ that some states of two interacting particles 
(TIP) in a random potential may be less strongly local- 
ized than each, particle separately has recently caused 
much interestEjo. This idea can be understoocL-&dthin 
the scaling theory based on the Thouless picturet^l. Ac- 
cording to this, the derealization of the particular TIP 
states from one block of size L\ (the one-particle lo- 
calization length) to the neighboring block follows from 
their having an interaction dependent Thouless number 
92 ~ E C 2/A2 much larger than that of the single-particle 
states g\. Here, A2 is the two-particle level separation at 
the given energy and the^corresponding Thouless energy 
E C 2 was first identifiedEj to the interaction dependent 
decay rate between neighboring blocks. As in the single- 
particle case, the multiple role played by the Thouless 
energy, as discussed above, immediately suggested that 
this TIP Thouless parameter will also be relevant to the 
level correlation problem. This is based on the general 
qualitative picture. While the TIP-spectrum without in- 
teraction contains only hidden one-particle correlations 
appearing on energy scales larger than Ai and is close to 
uncorrelated levels on lower energy scales, the interaction 
re-establishes the universal Wigner-Dyson rigidity up to 
the energy E C 2 = Ejj which depends on the strength U 
of the interaction. In the localized regime {L > Li), this 
was formally described by a nonlinear er-model for the 
TIP problem, as presented in Ref.E2l. The latter gives 
a theoretical foundation for the scaling picture for TIP 
on equal footing to that for non-interacting particles. In 
the metallic regime (L < Li), a study of the TIP-level 
statistical^ confirmed that Ejj gives also the characteris- 
tic energy scale up to which the TIP-spectrum exhibits 
the universal Wigner-Dysonjigidity. This was qualita- 
tively explained by mappings the TIP-Hamiltonian onto 
a Gaussian matrix model with preferential basis. 

We see that the scaling properties for interacting par- 
ticles can thus be studied via the spectral correlations of 
their levels in the metallic regime. This is an extremely 
useful observation. The metallic regime is easier to study 
both analytically, where reliable methods exist, and nu- 
merically. In the latter case, the necessity to go to very 
large system sizes larger than the localization length with 
weak disorder in low dimensions is eliminated thereby. 
Since the study of two interacting particles is only the 
first step towards the treatment of a more realistic many- 
body system, it is highly desirable to increase the number 
of particles. Even a modest program of going from two 
to three, four and larger numbers of interacting particles 
can be best accomplished by analyzing the Thouless pa- 
rameters in the metallic regime for rather small system 
sizes. This is the strategy we adopt in this paper. 

It was mentioned before (footnote 1) that when a state 
is coupled to a quasi-continuum, the golden rule expres- 
sion for its width starts to be valid only when the cou- 



pling is strong enough, or the density of the final states 
is high enough. The crossover between perturbative mix- 
ing (Rabi regime) and effective decay in fact occurs when 
the typical matrix clement of the coupling becomes Iff&Qf 
than the mean level spacing of the accessible statesu'LlliS. 
An equivalent condition is that the golden-rule width be 
larger than the final level spacing. This very general 
crossover, which becomes a phase transition in the appro- 
priate "thermodynamic limit" is the essence of dereal- 
ization in the usual Thouless scaling theory for ap-single 
particle. It applies to two-particle delocalizationEj and 
should likewise describe derealization for-JV particles. 
The Hilbert-space transition found in RefO is another 
example. In this case one gets a proper transition by 
the hierarchical coupling to higher and higher numbers 
of quasi-piarticle excitations. In the work presented here, 
as in Ref.E3, this transition is observed numerically as a 
function of the interaction strength U. When U is weak, 
it can couple only a few very close quasi-degenerate states 
and leads at most to Rabi-type oscillations between adja- 
cent levels. When U is larger than a certain threshold U c \ 
(or at larger excitation energy), many non- interacting 
states are coupled and Fermi's golden rule describes 
the spreading width of a non-interacting state over the 
(quasi-continuum) of other non-interacting states which 
are nearby in energy. U c \ is also the crossover interac- 
tion between Poisson and Wigner-Dyson statistics for the 
spectral fluctuations. We show that the same U c \ signi- 
fies also the breakdown of the perturbation theory in U . 
Above a higher threshold U C 2, the states are ergodically 
mixed and is suggested to increase like U N 

Most of this paper will be concerned with the three- 
particle problem. Simple analytical arguments will be 
presented for the behavior of the spectral correlations 
in a small diffusive quantum dot, and compared to a 
numerical study. Thus, we work in this paper in the 
metallic regime and do not directly study the derealiza- 
tion in real space for stronger disorder, when the one- 
particle states are localized. However, some remarks will 
be eventually made on the implication of this picture to 
interaction-induced derealization in real space and on 
its generalization to quasi-particle excitations in a de- 
generate metallic Fermi system. In particular, the basic 
derealization mechanism discussed in the original loca- 
tor expansion of Anderson was the divergence of per- 
turbation theory around the initial localized states. It 
will be shown that a seemingly analogous divergence can 
be identified in the perturbation theory in the interac- 
tion, around the noninteracting states. Here, this diver- 
gence signifies (as is also true in the Anderson localization 
case) the onset of Wigner-Dyson correlations in the full 
spectrum, where the basis of noninteracting eigenstates 
becomes well-mixed due to the interactions. A similar 
process appears in the recent workllll of Altshuler et al., 
using an analogous expansion for the quasiparticle excita- 
tions in a degenerate Fermi system, decaying by emitting 
electron-hole pairs. In the three cases of the Anderson 
derealization in real space and the derealization pro- 
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cesses found due to interactions in the Hilbert space of 
wavefunctions, the basic condition for the transition is 
very similar. It demands that the matrix element of a 
state to the "first generation" of states directly coupled 
to it by the interaction, be comparable to the level spac- 
ing of the latter states. 



sample size L. 7io, which contains the kinetic and poten- 
tial energy of the particles is combined with a two-body 
interaction of the form 



(i) 



ijkl 



II. JV-BODY HAMILTONIAN IN THE FOCK 
BASIS 

In the presence of interactions, it is convenient to con- 
sider the iV-body system in a certain Fock basis. Since 
we use this terminology in a slightly unusual way, let 
us make precise what we mean by Fock basis. We con- 
sider the one-particle states which take into account ex- 
actly the kinetic energy, the random electrostatic poten- 
tial seen by the electron, the chaotic or the integrable 
dynamics yielded by the boundaries in a ballistic billiard 
etc., and we use the exact one-particle states to build up 
the Slater determinants which we refer to as the Fock 
states. Therefore, by Fock basis we just mean the eigen- 
basis of the iV-body Hilbert space in which the system 
Hamiltonian is diagonal at U = 0. 

Moreover, we do not focus on the low excitation ener- 
gies, (i.e. on the restricted space available from the Fermi 
vacuum by successive applications of quasi-particle cre- 
ation operators) but rather to higher energies in the bulk 
of the TV-body spectrum. Therefore, in contrast to Ref.EJ, 
the parameter in our study is not the excitation energy 
of an extra quasi-particle above the Fermi sea, but the 
strength U of the interaction, at a given total energy 
chosen close to the band center of the iV-body spectrum. 
Another important difference between this study and the 
problem considered in Ref.tll is that we have in mind the 
dilute limit [], where the number of "real" particles is ar- 
bitrary, but nevertheless of zero density. Therefore, we 
have not in this study a Fermi vacuum from which an ar- 
bitrary large number of quasi-particles can be indefinitely 
created. 

In this Fock basis, the Hamiltonian with interaction 
is a random matrix with preferential basis. For the 
sake of simplicity, we assume that the one-particle states 
are more or less uniformly extended inside the sample 
(no one-particle localization); i.e. the Hamiltonian with- 
out interaction Tio is a sum of one-particle Hamiltonians 
which can be described by a random matrix being sta- 
tistically invariant under the orthogonal transformation 
0(M). M = L d is the number of considered sites for a 



t This dilute limit strictly means that N/M — ► when 
M — > oo, where M = L d is the number of sites in a tight- 
binding model. This should be distinguished from the finite 
density limit where N/M is a constant and the Fermi system 
is degenerate at zero temperature 



where |i) denotes the Wannier function localized at the 
site, Cj /cf destroys/creates a particle on site i, and 



;th 



U = 



E 



ij)uij{ij\ 



(2) 



is a local interaction of strength measured by a param- 
eter U. We will give estimates assuming on-site inter- 
action in the following, but the interaction should be 
obviously extended to nearest neighbors in the case of 
spinless Fermions. We write the Hamiltonian in the 
basis of the n tot = M\/(N\{M - N)\) Slater determi- 
nants (Fock states), which are antisymmetrized products 
\A n >= \a±, . . . , a 7v > (n = 1, . . . , ntot) of one-particle 
eigenf unctions \aj >. 7io is a diagonal matrix with the 
different possible sums Yl7=i €a i °f one-particle energies 
as entries. The interaction term H. ln t yields a matrix with 
entries 



<a[ . . . a' N \H- m t\ai ■ ■ - a N >=E IT ^ a 'i a i\ Q 



IJ \i^IJ 



(3) 



where 



M 

E 

p,p'=i 



(4) 

tjj ai (p) denoting the amplitude of the wave-function in 
the one-particle eigenstate \aj > at the site p. 

The existence of the interaction yields two effects that 
we consider separately. The diagonal matrix elements 
of Hint shift the location of the iV-body levels, an ef- 
fect which is predominant for small system size and 
large strength of the interaction, and which can yield 
an i mporta nt re-arrangement of the spectrum (see sec- 
This situation is shortly described in the 
ftion, and has been extensively discussed 
in the case of Anderson insulators with 
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following subs 
by Kamimurac-S 
a very small localization length. The off-diagonal matrix 
elements of H.- m t give rise to hopping among certain Fock 
states, and thus to derealization in the Fock basis. In 
this study, we mainly focus on the description of this in- 
teraction induced derealization in the Fock basis, in the 
limit where the second effect dominates the first. This 
derealization in the Fock basis is a generic effect of the 
interaction which should not be confused with derealiza- 
tion in real space. It is only when the one-particle states 
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are themselves localized in real space that Hilbert space 
delocalization may-result in delocalization in real space. 
A recent analysis[jp3 of the sensitivity of the energy lev- 
els to a change of boundary conditions has stressed this 
difference. 



A. Diagonal matrix elements of the interaction and 
large [/-limit 

For very large U and small system size L, the previ- 
ously defined Fock basis is no longer appropriate. It is 
more instructive to consider the Fock basis built of the 
on-site orbitals, and not of the one-particle eigenstates. 
The kinetic energy, and not the interaction, can then 
be treated perturbatively. In this basis the TV-particle 
states without kinetic terms can be classified according 
to the number of next-neighbor configurations, for a next- 
neighbor interacti on. Th is limit will be discussed in more 
detail in section VIII D| , where numerical results show 
that at U ~ 15, the spectrum of three spinless Fermions is 
split into three separated bands, with a density of states 
approximately given by the sum of three Gaussians cen- 
tered at E = 0, U, and 2U. The weights of those three 
bands are directly related to the number of next-neighbor 
configurations in the on-site Fock basis states. 

However, for U ~ 1 we are far from seeing interac- 
tion induced gaps in the spectrum. We then assume 
that the diagonal matrix elements of the Hamiltonian in 
the Fock basis built of the one-particle eigenstates, are 
mainly dominated by the one-particle contributions, and 
that the effects coming from the interaction can be ne- 
glected for those elements. Therefore, we consider only 
the delocalization in this basis, which results from the 
off-diagonal terms. 



B. Off-diagonal matrix elements of the interaction 

For TV > 3, one can see from Eq. (||) that there are 
only non-zero matrix elements between Fock states hav- 
ing TV — 2 one-particle states in common. This means 
that the Hamiltonian (|l]) is a sparse matrix in this Fock 



basis, with a strongly preferential basis_T-his property 
was not present in the former studiesoliSEJ for TV = 2, 
and a straightforward generalization of the former re- 
sults to TV particles would require TV-body interactions. 
The two-body form of the interaction introduces specific 
problems for TV > 3, which has been recently discussed in 
Ref.cJ'EZl. Moreover, there is a large degeneracy of these 
non-zero terms. For instance, when TV — 3, all the ele- 
ments < a 1 a.2(3\'Hi n t\a' 1 a' 2 (3' >= Q ai a 2 a' ia 'Jpp' are the 
same for all of the one-particle states j3 = j3' '. 

The form of the distribution of the degenerate non-zero 
off-diagonal terms (Eq.[|) is by itself a non-trivial one- 
particle problem. If the underlying classical one-particle 
dynamics is diffusive, as in a disordered metal, ballistic 
chaotic, as-,p, a billiard, or integrable, one gets different 
estimatestlrea for the magnitude of the interaction matrix 
elements. For simplicity we will use the very rough ap- 
proximation of uncorrelated one-particle wave-functions, 
with amplitude of the order l/VM on each site with a 
random sign. This corresponds to a one-particle Hamilto- 
nian being statistically invariant under orthogonal trans- 
formations (O(M) invariance assumed in standard. Ran- 
dom Matrix Theory). As pointed out in Ref.liilS, this 
evaluation of the interaction matrix elements only repro- 
duces^] the zero wave mode contribution of a diffusion 
process. This gives Q a ' 1 a' 2 a 1 a 2 ~ ±<5t yP ~ ±U/M 3 ^ 2 
for on-site interactions. We use this approximation for 
the simplified theoretical picture that we present before 
the numerical study. This is because we want to compare 
our predictions to simulations on disordered systems with 
too small sizes to have one-particle diffusion. The more 
detailed description of the off-diagonal matrix element 
will modify the quantitative dependence as a function of 
the system parameters, but will not change the general 
scheme of the effect of the interaction on the spectral 
correlations. 

Moreover, we will neglect the energy dependence of the 
TV-particle density p^, taking p N x — Ajv ~ B\/M N . An 
denotes the typical TV-particle level spacing in the bulk of 
the TV-particle spectrum and B± is the bandwidth (one- 
particle kinetic energy scale). 



*In RefM it is also shown that the relation between the 
spectral rigidity and the level curvature is not direct for the 
TV body problem. For TV = 1, the original definition of the 
Thouless number, in terms of the curvature of the energy lev- 
els, coincides with our definition, based on the spectral rigid- 
ity. For TV = 2, these two definitions are not the same in the 
metallic regime, but the (spectral) Thouless number can be 
expressed in terms of another (topological) curvature, assum- 
ing that distinct Aharonov-bohm fluxes can be associated to 
the distinct particles. For TV > 3, the study of those different 
curvatures is postponed to a further study. 



^We note two important modifications to our rough approx- 
imation for the matrix elements. For energy transfer smaller 
than the Thouless energy E c , its magnitude is enhanced to 
the order of ±Ai/<?i when g\ > 1. In one dimension (gi < 1), 
it was pointed out by Ponomarev and SilvestrovCJ that there 
are important modifications of the-4natrix elements. More 
precisely, as recently shown in Ref.tJ, when the one particle 
states are localized, the fluctuations of the interaction ma- 
trix elements are so large that the effective density of directly 
coupled Fock states becomes multifractal. 
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III. SPREAD WIDTH OVER STATES DIRECTLY 
COUPLED BY THE INTERACTION 

For TV = 2, the full Hamiltonian can be modeledEl by 
a Gaussianj-ensemble of random matrices with a prefer- 
ential basisEil. The structure of the projections Cp n =< 
'fyplAn > of the many-body eigenfunctions \*f?p > (la- 
beled by (3 = 1, . . . ,n to t) onto the Fock,-steies \A n > is 
well described by the Breit-Wigner formllaEJ 



{\C(>n\ 



■2n[(E -E n y 



rl/4] 



(5) 



where the brackets denote ensemble averaging, and the 
spread width 



U 2 M 2 

T 2 =2,{Q 2 ) P2 *2,— — 



(0) 



increases with the interaction according to Fermi's golden 
rule. This means that (for T 2 > A 2 ) an eigenfunction 
\9p >= J2n Cf3n\A n > has significant projections on typ- 
ically T2/A2 Fock states. 

For TV > 3, this can be generalized to a spreading 
width Tff oc (Q 2 }pff where pjy is the density of the 
iV-body Fock states directly coupled by the interaction. 
For spinless Fermions, one has n to t = M\/(N\(M — TV)!) 
Fock states and the number of Fock states directly cou- 
pled by the two-body interaction is = N(M — N) + 
TV(TV - 1)(M - N)(M - TV - l)/4. In the dilute limit 



oc M 



f2 



TV < M, one finds n tot oc M N /Nl while 
is much smaller than n to t- Assuming uniform densities, 
this means that the effective level spacing between 
Fock states directly coupled by the interaction is of or- 
der AS* Bi/rigg . For very few particles, AS* « A2, 
but we emphasize that this approximate relation becomes 
uncorrect if TV is large, mainly in the finite density case 
where N oc M. 

A quantity closely related to the local density of states 
is the participation ratio R — (X)"=i |C/3n| 4 ) -1 , which 
gives the number of Fock states mixed by the interac- 
tion. Using the structure @ of the eigenfunctions at 

T 2 > A2, one can get the estimate R ~ tvT^ /A^ « 
2n 2 U 2 (M 3 B 2 )~ 1 (n ( c d s !) 2 . Therefore we expect to find 
R oc U 2 , since the contribution of the states directly cou- 
pled by the interaction will dominate for small U . 

Therefore, the first observable effect of the interaction 
will be the broadening of a Fock state over r^/A^ 
other Fock states separated by a characteristic scale 



A 



(d) 
N 



Bi/n 



(d) 



1 eff 



This spreading width is pro- 



portional to U 2 , but does not characterize the coupling 
of the original Fock State to the TV-body spectrum. In 
this spread width, there are many Fock states (of a den- 
sity ps = I/A3 for TV — 3) which are not directly coupled 
to the original Fock state at this order in U. This is the 
major difference between N = 2 and TV > 3: For TV = 2, 



the width T2 which characterizes the local density of in- 
teracting states in the Fock basis is directly related to the 
spectral statistics: the energy scale E\j up to which the 
spectpum exhibits the universal Wigner-Dyson rigidity is 
giventS by this spread width T2 , provided T2 > A2 . For 
N > 3, even when rff > An, the levels separated by 
Ajv are not necessarily coupled and can be statistically 
independent. For the level repulsion at the scale Ajv, the 
spreading width Tff does not provide a relevant energy 
scale. The A^-particle Thouless number gjq is not given 

by r^/A^, when TV > 3. 

In the following, we discuss what should provide this 
relevant energy scale for the spectral statistics of con- 
secutive TV-body levels, with a density P m = 1/Ajy, and 
thus the relevant gM- We consider first the case TV = 3. 
The generalization to an arbitrary number TV of particles 
is straightforward, as far as we are in the dilute limit. 



IV. PERTURB ATIVE REGIME (U < U C i) 

The spectrum without interaction contains only one- 
particle correlations appearing on energy scales larger 
than Ai and is close to uncorrelated levels-pn lower en- 
ergy scales. The interaction re-establishesEj the univer- 
sal Wigner-Dyson rigidity up to an energy Ejj which 
depends on the strength of the interaction. When U 
is weak, it can couple only quasi-degenerate states and 
leads at most to Rabi-oscillations between adjacent lev- 
els. 

For TV = 2, E v oc \U\, while one finds E v = T 2 oc U 2 
at larger U when g 2 = T2/A2 > 1, i.e. when many Fock 
states are coupled rajjid Fermi's golden rule applies. This 
can be understoodla from the following arguments. For 
very weak U (T 2 < A 2 ) only the coupling between two 
Fock states with a separation < A2 is relevant. This re- 
stricts the problem to the analysis of a solvable 2x2 ran- 
dom matrix, with diagonal terms typically much larger 
than the off-diagonal coupling term. A model with in- 
dependent Gaussian entries, the variance of the diagonal 
entries being much larger than_these of the off-diagonal 
term, was exactly solved in Ref.El Since a 2 x 2 real sym- 
metric matrix can be diagonalized by a rotation of angle 
9, it is easy to write the probability distribution in terms 
of the two eigenvalues and of 8. Integrating over 9 gives 
the joint probability distribution of the two eigenvalues. 
One finds for this 2x2 matrix model that E\j is given 
by the absolute value (r.m.s) of the off-diagonal term. 
It was interpreted in terms of Rabi oscillations between 
two Fock states at typically A2 away from each other 
in energy and the range E\j of the level repulsion was 
identified)**] with this Rabi frequency. For energy separa- 



**A related effect is for the one-particle prob- 

lem, where the Thouless sensitivity to boundary conditions is 
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tion e < Ejj , the consecutive levels repel each other as in 
standard random matrix theory, while their fluctuations 
are uncorrelated for e > Ejj. 

For N = 3, we denote by \A n >= |aia 2 a3 > the three- 



particle Fock states, of energy E„ = e c 



, , and 



we consider two energetically nearby Fock states | A n > 
and | A n / >: i.e. with E n — E n i ss A3, the three-particle 
level spacing. For a weak interaction, as pointed out 
by Shepelyansky and SushkovcJ, the effective matrix el- 
ement Uf of the interaction between those two consec- 
utive Fock states can be estimated using perturbation 
theory. It is only in second order that one gets a non- 
zero contribution resulting from terms like 



E 



< A n \U 12 \A n „ >< A n „\U 23 \A n , > 



E„ 



E„n 



(7) 



where particle 1 interacts with particle 2, then particle 
2 with particle 3. Since we have a two-body interac- 
tion (< A n \U 12 \A'^ >=< a 1 a 2 \U 1 2\a'{a'^ > £ a3 «»), the 
summation over n" is reduced to a sum over the single- 
particle quantum number a 2 . This sum is of the order 
of its largest term, i.e. of a term with an energy denomi- 
nator of order Ai, the one-particle level spacing, and not 
A3. This eventually gives for the effective matrix clement 
which couples two consecutive Fock states a magnitude 
of order 



U: 



off 



±- 



typ 



Ai 



± 



U 2 



A/ 3 Ai 



(8) 



if one takes for U typ our simple estimate 
Q typ = ±U/M 3 / 2 . 

Therefore, in this perturbative regime, Fock states at 
from each other are coupled by a matrix 



A( d ) 

an energy A 3 



element Uf' « U typ ~ ±U/M 3 > 2 while Fock states at A 3 
from each other are only coupled by Uf cx ±U 2 yp /Ai w 
±U 2 / (M 3 Ai). From this, we draw two main conclusions 
for the three particle problem that we extend to an arbi- 
trary number N of particles. 



A. Hierarchical structure in the Fock basis 
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FIG. 1. Structure of the eigenfunction at U < U c i in the 
Fock basis 

broadened over a density pz = A3 , with a Breit-Wigner 



shape characterized by a width 
perturbative regime, r 3 i} < A 3 when U < U cl 



(Uf) 2 A3 1 . In the 



B. Rabi frequency and Wigner-Dyson rigidity 

When U is so small that the broadening is smaller 
than A3, the levels are essentially uncorrelated. How- 
ever, level repulsion occurs for energy scales smaller than 
Ejj w |£/f ff | < A3. The reason for this is a straightfor- 
ward generalization of the argument given for N = 2: 
The energy level correlations come from the very small 
coupling terms between Fock states which are nearest 
neighbors in energy. This reduces the complicated prob- 
lem of a very large random matrix with a few non-zero 
small off-diagonal terms to the solvable problem of a 2 x 2 
matrix with an effective off-diagonal term of magnitude 
U§ . Rabi-oscillations between the two coupled diagonal 
terms occur. Their frequency, of order |J7f ff |, character- 
izes also the scale Ejj below which the universal level 
repulsion occurs. For TV = 2, this gives Ejj oc \U\ (direct 
coupling), while E v w \Uf \ cx |J7 2 | for N = 3. 

Let us consider now the case where N = 4. As for 
the three particle case, it is sufficient to go to the sec- 
ond order in perturbation theory for having a non zero 
effective matrix element coupling two consecutive Fock 
states. One has: 



The states in the Fock basis have a very particu- 
lar hierarchical structure, as sketched in figure U A 

~ (d) - 1/A$ d) « 



Fock state is broadened over a density p\ 



l/Aj of neighboring Fock states. This broadening has 
a Breit-Wigner form characterized by a width s» 
C/ 2 /(^ 3 A 3 d) ). The projections over the neighboring 
Fock states at A 3 d ' ) away from each other are themselves 



proportional to the square root of the Landauer conductance 
when gi < 1 and to the Landauer conductance when gi > 1. 



m 



E 



< A n \U 12 \A n „ >< A n „\U 34 \A n , > 
E n — E n » 



(9) 



where the \A n > are now the Fock states for 4 particles. 
The difference with the case where N = 3 is that the 
summation over n" is now totally suppressed. This yields 
a smallest possible denominator E n ~E n i> of order Bi and 
not Ai. One finds J7| ff w ±U 2 yp / Bi. Similarly, one finds 
Uf w ±U* yp /(B 1 A 1 ), Uf w ±U? yp /(B 2 )... and the 
general expression is given by Uf w ±Bi(U typ / Bi) p for 
an even number N of particles with P = N/2, and Uf w 
±({7 typ / J Bi) p ( J B 2 /Ai) for an odd number N with P = 



G 



(N + 1) /2. This gives us the energy range Ejj pa |Z7|p | of 
the level repulsion for weak interaction and arbitrary N. 



V. CROSSOVER FROM POISSON TO WIGNER: 

U = Uci. 

When r<° re (Uf) 2 /A 3 Pa A 3 , the Wigner-Dyson 
rigidity is established on the scale A3 for the three parti- 
cle case. This defines the first interaction threshold U c \ 
where a sharp cross-over from Poisson to Wigner should 
be observed in the distribution P(S) of the spacings be- 
tween consecutive energy levels. For an arbitrary number 
N of particles, the first threshold U c \ is given by the con- 
dition 



9n 



Ml 



1 . 



(10) 



Using the estimate for Ufj given by the first non-zero 
order in the perturbative expansion in U , and assuming 
that N is small enough to have A 2 pa Af* , (i.e. neglecting 
a factor of order 1/N 2 ), one finds the general relation 



9N 



fill 



where P = N/2 when N is even and P = (N + l)/2 
when N is odd. Therefore, the interaction U c i where the 
two particle Thouless number g 2 , given by (U/Bi) 2 M in 
our estimate, is of order one (i.e. Ut yp ~ A| ff ), does not 
signify a Poisson- Wigner cross-over for the two particle 
case only, but also for the iV-body spectrum, as far as 
we stay in the dilute limit. We note that our general 
relational implies g 3 ~ g%, in agreement with Ref.E3. 

U c i can also be understood from the parametric mo- 
tion of the energy levels when U increases. Let us con- 
sider again the case where N = 3. Level repulsion starts 
to be efficient when a certain relative characteristic en- 

("3") 

ergy change (U) of the levels due to the interaction 
is of the order of the mean level spacing. The displace- 
ment resulting from the diagonal matrix elements of the 
interaction induces a parallel motion of the levels and is 
therefore irrelevant. For the real part of the self-energy 
Sr(C7) of a given Fock state, perturbation theory gives 
different terms involving loops starting from the consid- 
ered Fock State and visiting one, two, three and more 
Fock states being directly coupled among each other by 
the two-body interaction before returning to the start- 
ing point. The term of order U 2 (loop visiting a single 



"In Eq. the Thouless numbers are defined by the en- 

ergy scale below which level repulsion occurs, in units of the 
mean level spacing Ajv. This definition, valid in the pertipi 
bative regime only, differs from the definition used in RefH3 
(g N = F^/A N w (t/^f /Ajv) 2 ). However, this subtlety (see 
the previous footnote) does not matter for Eq. (111]). 



Fock state) cannot be relevant since it involves an energy 
denominator of order A| , and not A3. The first term 
involving loops visiting Fock states at A 3 away from the 
unperturbed level is only given by the term of order U 4 
(loop starting from the Fock state and visiting three dif- 
ferent Fock states before return). For the real part, this 
gives terms like 



< A n \U l2 \A n > >< An>\U 23 \A n » > 

n'n"n'" 



(E n — E n >) 

< A n „\U 23 \A n m >< A n n,\U 12 \A n > 



{E n -E n „){E n ~E n ,„) 



(12) 



These sums are of the order of the term having the small- 
est possible denominator, e. g. [E n — E n i>) pa A 3 , which 
fixes the state | A n ii > and thus suppresses the sum- 
mation over n" . Then, the two-body character of the 
interaction yields for the energy differences E n — E n > and 

E n — E n "i , with E n and E n " fixed, a smallest possible 

( 1) 

value of order Ai. One can see that ££'(£0 ps A3 pre- 
cisely when 53 = I3 /A3 pa 1, for U = U c \. It is straight- 
forward to check that it is also at U = U c \ that the term 
of order U A of the perturbative expansion of the imagi- 
nary part of the self-energy is of order A 3 . Moreover, this 
argument can be easily extended to an arbitrary number 
iV of particles, with a conclusion in agreement with these 
previously presented in order to obtain U c \. 



VI. BREAKDOWN OF THE PERTURBATION 
THEORY AND NON-ERGODIC 
WIGNER-DYSON REGIME: Uci < U < U C 2 

The argument that we propose to characterize the rele- 
vant Theailsss numbers g n is reminiscent of a locator ex- 
pansioiS'E^I "a la Anderson" , where the self-energy of the 
Fock states in the presence of the hopping terms yielded 
by the interaction is evaluated using perturbation theory. 
For one-particle localization, the breakdown of this per- 
turbative expansion in the basis built of the site orbitals 
was related to a metal-insulator transition due to dereal- 
ization in real space. One knows from the scaling theory 
too that this transition occurs for g\ sa 1 . This is a closely 
related consideration which has led the authors of Ref.U 
to conjecture that interactions should give an Anderson 
transition in Fock space, for a critical value e* of the ex- 
citation energy of the extra quasi-particle injected above 
the Fermi sea. This led us to determine up to what maxi- 
mum value of U the relevant Thouless numbers gjsr can be 
given by perturbation theory. If one follows Shepelyan- 
sky and SushkovEH, who assume that perturbation theory 
remains valid when 53 > 1 for N — 3, the relation g 3 pa g\ 

gives (73 = r 3 ^/A3 cx C/ 4 above U c \. However, in our nu- 



merical study (see section VIII) we observe above U c \ a 
|J7|-increase of g 3 , following the perturbative i7 2 -increase 
(Rabi regime). Moreover, all the quantities calculated for 



7 



three spinless Fermions (local density of states, partici- 
pation ratio, spectral rigidity) have not given any trace 
of a [/ 4 -proportional behavior when U > U c \. 

This leads us to suspect that perturbation theory can 
only be used up to U c \ 1 where level repulsion is estab- 
lished at the scale A3. Above this threshold, the rele- 
vant T3 does not coincide with the perturbative estimate 



This breakdown of the perturbation theory above 



U c i, when Ut yp > A!j ff , can be shown if one evaluates 
the effective matrix element coupling nearby Fock states 
at higher orders in U . Let us present the argument for 
N = 4. The coupling term of order U 2 was found to be 
[7f ff « ±/7 t 2 yp /i?i. This corresponds to a process where 
particle 1 interacts with particle 2, then particle 3 with 
particle 4. A term of order U 3 is given for instance if 
particle 3 interacts with particle 4 once more. 



Uf (order 3) « 
<A n \U 12 \A Ti 



>< .An" I #34 1 A."' >< A. 



(13) 

\U 3i \An>> 



(E n — E n n)(E n — E n m) 



The smallest energy denominator is the product of two 
energy differences. The hrst one is of order B±, but one 
has extra degrees of freedom for the choice of two one 
particle quantum numbers of \A n m > since particle 3 in- 
teracts with particle 4 two times. This reduces the second 
energy difference to a smallest value of order Af , and 
gives Z7f ff (order 3) w Uf (U typ /Af). This means that 
the terms in U 2 and in U 3 of the perturbative expansion 
are of the same order for Ut yp ~ A| . The generalization 
is straigthforward: For arbitary N, one finds 



Uf (order p) ra Uf 



'N 



U 



typ 

v cff 



p—p 



(14) 



for the terms of order p > P, where P is the order giving 
Uf . This indicates that the sum to all orders in U does 
not converge above U c ±. Similar considerations, used here 
to evaluate to all orders the effective matrix element cou- 
pling Fock states which are nearby in energy, can be given 
for the self-energy. From this emerges the general result 
that the perturbative sum in U has terms with similar 
magnitude to all orders when U typ ~ A| e . One can then 
conclude that perturbation theory breaks down at the 
Poisson-Wigner cross-over in the ./V-body spectrum. If 
one compares our conclusions with the ones presented in 
Ref.E3, we also find that the threshold U c \ corresponds 
to gi ~ 1 for very few particles. However, in contrast to 
the conclusions of RefEO, we have shown that the rela- 
tion f/3 sa g 2 for N — 3 cannot be extended for g 2 > 1 
and holds only when g 2 < 1. For finite N in a finite size 
system, the spectral statistics exhibits a cross-over from 
Poisson to Wigner. When N increases, this cross-over 
should become sharper and sharper to eventually give 
a real transition. It is natural that such a transition is 
accompanied by a breakdown of perturbation theory of 
the self-energy of the Fock states, for /7 typ ~ A§ ff , as it 
happens in the locator expansion of the self-energy for 




FIG. 2. Structure of the eigenfunction at U c i < U < U C 2 
in the Fock basis 

the one-particle problem at g% « 1. In both cases the 
condition is that the matrix element of a state to the 
"first generation" of states directly coupled to it by the 
interaction, becomes comparable to the level spacing of 
the latter states. A similar picture and derealization 
condition applies also to the analogous transition in the 
Hilbert space of different numbers of excited ouasiparti- 
cles in a degenerate Fermi system, suggestedEJ recently 
by Altshuler et al.. 

As mentioned above, g 3 does not increase as U A when 
33 > 1 (see figure ||). We are however not able to explain 
the observed linear increase for U > U c \. Nevertheless, 
one can say that the interacting states for U > U c \ should 
not be ergodic in the energy window where they are 
broadened, but still have a structure as sketched in figure 
|. When T3 = 53 A3 is much smaller than A^ « A| ff , 
there are still many Fock states inside the energy width 
where an interacting state has essentially a zero pro- 
jection. This will disappear only at a second threshold 
U C 2, characterized by the condition: T3(U C 2) = A 3 d ^. 



VII. ERGODIC WIGNER-DYSON REGIME AND 
SELF-CONSISTENT THEORY: U > U C 2 

We present here a conjecture for the regime of rather 
large interactions, when a Fock state is well coupled by 
the interaction to the three-body spectrum. We consider 
again the case where N — 3 and a sufficiently strong in- 
teraction for having T3 > A^, but nevertheless small 
enough for not being in the large [/-limit dominated by 
the diagonal terms of the interaction. For U > U c %, one 
can assume that the interacting states are not unambigu- 
ously related to the previous hierarchy of Fock states, 
but are closer to random mixtures of 33 consecutive Fock 
states, each of them contributing with a projection of 
random sign and of typical amplitude of order l/^/g^. 
In other words, one has a simpler case where all the 
Fock states in an energy window F 3 are now well cou- 
pled, and remain decoupled from the other Fock states 
outside this window. This is what we mean by "ergodic 
Wigner-Dyson regime" , where the sparsity of the 3-body 
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FIG. 3. Structure of the eigenfunction for U > U C 2 in the 
Fock basis 

Hamiltonian yielded by the pairwise character of the in- 
teraction becomes essentially irrelevant. A self-consistent 
evaluation of c/ 3 becomes possible if one assumes that an 
interacting state \^>p > has the following structure in the 
Fock basis: 



i^>= E 



(15) 



For N — 3, we calculate the interaction matrix element 
Qp'P between two states \^fp > and \iffpi > at nearby 
energies, each of them being superpositions (with am- 
plitude Cf lQ2Q3 « l/Vsi) 01 53 = r 3 /A 3 Fock states 
(A 3 = Sx/ntot) for L < L x : 

Q/3'P = E ^ojai^oioaoa < a 'l a 2«3 l^int > 

a l a 2 a 3 
Q l a 2 a 3 

(16) 

where the sums have to run over g 3 basis states. This 
means, that each of the a-summations runs over g^ 3 
values. As can be seen from (^), the matrix elements 
between the Fock states contain three terms of the form 
Qa' I a.' J a.ia.j&a.' K a,K wnere {7, J, K} are the different cyclic 
permutations of {1,2,3}. The Kronecker-(5 reduces the 
relevant summations occurring in ([l6]) to 5 summations, 

1 /3 

each of them running over roughly g 3 values. Thus, 

5/3 

Qp>p consists of a sum of 3<? 3 terms of typical size 
Qtyp/ff3 with random sign, which yields the result 



(\Qp>p\ 2 ) « 3 5: 



5/3 Qty p 
3 „2 



W 2 



9s M 3 gl /3 



(17) 



Plugging this into a Fermi golden rule evaluation of the 
spread width 



r 3 = 



M\Qp>p?) 



one obtains 



53 



T 3 (6tt) 3 / 4 / u 



A 3 ~ M 9 / 4 V nt0t ^i 



3/2 



(18) 



(19) 



At U > U C 2, we therefore expect the decay width of the 
eigenf unctions and the participation ratio R to increase 
like U 3 / 2 . 

For N spinless Fermions, this ergodic Wigner Dyson 
regime is characterized by as it can 

be seen from a straightforward generalization of the self- 
consistent argument presented above for N = 3. 



VIII. NUMERICAL STUDY OF THREE 
SPINLESS FERMIONS 

In order to illustrate our theory, we have performed 
a numerical study of the many-body eigenstates and 
eigenenergies for three spinless Fermions in a disordered 
cubic lattice. 



A. Numerical model and characteristic scales 

For the numerical simulations, we use a three dimen- 
sional tight binding model on a cubic lattice containing 
3x3x3 sites. The disorder and the hopping terms are de- 
scribed by the usual Anderson Hamiltonian with on-site 
potentials drawn from a rectangular distribution of width 
2W with W — 2 and nearest neighbor hopping terms 
V = 1 which set the energy scale. We use rigid boundary 
conditions in all three directions. In addition, we use a 
two-body interaction of the form (|l|) with Ujj = U when 
the sites i and j are nearest neighbors on the lattice and 
Uij = otherwise. 

In such a cube, there are M = 27 one-particle states 
with a typical density of p\ = 1/Ai w 4 in the cen- 
ter of the band. For spinless Fermions, this leads to 
M 2 = M (M - l)/2 = 351 two-particle states and 
n tot = M 3 = M(M - 1)(M - 2)/3! = 2925 three- 
particle states. For U ps 1 and in the center of the 
band, the density of the three-particle levels is about 
P3 = I/A3 w 270, while the density of two-particle 
levels amounts to P2 = I/A2 ~ 40. The density of 
three-particle Slater determinants directly coupled to a 
given state by the interaction is larger: with the number 

nfj = N(M-N)+N(N-l)(M-N)(M-N-l)/4= 900 
of non-zero off-diagonal interaction matrix elements in 



a line of H- m t, one finds p\ 



(d) 



1/A 



3 



1/A 



off 
2 



83. 



«iff/( n totA 3 ) 

For the analysis of the numerical results, we will 
slightly improve our estimates of the interaction matrix 
elements, taking into account that the interaction is not 
strictly on-site, but of range 1, since a particle on a given 
site can interact with another one when the latter is on 
one of the adjacent sites. The matrix element coupling 
two states of the Fock basis ([J) contains a double sum 
over the sites p and p' of the lattice. There are non-zero 
contributions to the sum whenever site p is a next neigh- 
bor of site p on the lattice. In the cube we consider, there 
is one site which has 6 next neighbors (NN), 6 sites with 
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5 NN, 12 sites with 4 NN and 8 sites with 3 NN. The 
mean number of next neighbors is Z — 4 and the sums 
in (0) run over a total number of ZM terms. Assum- 
ing the statistical invariance of the one-particle Hamilto- 
nian under orthogonal transformations, this yields a typ- 
ical size of the off-diagonal interaction matrix elements 
Qtyp ps ±UVZ/M 3 / 2 = ±0.01417. 

In the same way, one finds Q ai a 2 a 1 a 2 ~ 3ZU/M w 
0.45(7 for the diagonal terms of the interaction, which 
lead to a shift of the diagonal elements of the Hamilto- 
nian and conserve the sign of the interaction. The factor 
of three is due to the combinatorial factor which counts 
the number of different pairs out of three particles. For 
not too large U we neglect them as compared to the fluc- 
tuations of the diagonal elements (which are of the order 
of the band-width B 1 w 10) for U = 0. 

This must be carefully taken into account when es- 
timating statistical properties of these matrices. In the 
Fermi golden rule formula for the spread width of the lev- 
els (||) one must introduce the effective level spacing Ag" 3 ' 
of the directly available levels and finds the expression 



r<j d) = 2ir(Q 2 )/A 



off 
2 



2irU 2 ZM- 3 p ( 3 d) 



O.llt/ 2 . (20) 



The effective interaction between basis states which are 
coupled by second order processes only, can therefore be 
estimated to be 



u: 



eff 



± 



3\ZU 2 
A7 3 Ai 



leading to the spread width 



,(0 



oJ U f? 
Ztt — 



0.002CT 



0.0067<7 4 



(21) 



(22) 



in the perturbative regime where r 3 < A3. From these 
estimates, we get the first threshold U& ~ 0.85 for which 



typ 



AGff 

^2 ■ 



B. Structure of the wave functions 

We first concentrate on the structure of the eigenstates 
in the Fock basis. Examples are shown in Fig. [|. In each 
of the pictures, only one eigenstate \^Bp > at an energy 
Eg ~ is shown. Each point represents the overlap 
\C/3n\ 2 = I < */3|^n > | 2 with a Fock state \A n > and 
is plotted as a function of the energy difference Ep — E n 
between the eigenstate and the Fock state. 

It can be seen that in the case of two particles, al- 
most all of the Fock states which are in a certain energy 
range around the energy of the unperturbed eigenstate 
have a non-negligible overlap with it. For three particles, 
however, many very small values of the projections onto 
Fock states occur, even when they are quite close in en- 
ergy. However, it is difficult to observe the hierarchical 
structure of the three-particle states because the scales 
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FIG. 4. The projections |C/3„| 2 of a typical eigenstate 
> (Ep m 0) at interaction strength U = 0.5 (left) and 
U — 2.0 (right) on the basis of Fock states \A„ >. The data 
are for N = 2 (upper) and N = 3 (lower) spinless Fermions 
in a cube of 3 x 3 x 3 sites with interaction U between next 
neighbors and on-site disorder W = 2. 

A3 and Aj differ by a factor of 3 only in our case and 
because of the statistical fluctuations. 

From these overlap matrix elements, taking into ac- 
count several different realizations, we have computed 
the local density of interacting states in the Fock ba- 
sis (Wigner strength function). In spite of the fact that 
there are many very small values in the individual over- 
lap matrix elements, the average is well described by the 
Breit- Wigner form (||). Its spread width T3 is shown in 
Fig. U (lower data points) as a function of U. Therefore 
by T 3 here, we mean the total spread width extracted 
from the average local density of interacting states in 
the Fock basis, and not the partial spread widths 
and F 3 ^ introduced previously. T3 behaves quadratically 
down to rather low interaction values while T3 oc U 3 ^ 2 
above U ~ 1.5. 

For weak interaction, we are in the regime where the 
hierarchical structure of the eigenstates should be im- 
portant. The spread width is then dominated by the 

spread width Our estimates presented above give 

r 3 d) w 0.1 IU 2 which is the correct order of magnitude. 

From our theoretical considerations, we expect to ob- 
tain a regime in which the wave- functions are ergodic 
and the sparseness of the Hamiltonian irrelevant when 
U > U C 2 ~ 1.8. Taking into account the refined esti- 
mates of this section, the spread width is expected to be 
given by the self-consistent expression 



r - a ( 67rZu2 



3/4 



0.25[/ 3/2 . 



(23) 
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FIG. 5. Upper curve: Participation ratio 7? for 3 spinless 
Fermions in a 3 x 3 x 3 disordered lattice (V = 1, W = 2) with 
rigid boundary conditions and nearest neighbor interaction U. 
The dashed, dotted and solid line represent power law fits in 
the different regimes yielding R = 1617, 60U 2 , and 67U 3/2 , 
respectively. Lower curve: Spread width T3 characterizing 
the local density of states of 3 spinless interacting Fermions 
in the Fock basis. The dotted and solid line represent power 
law fits in the different regimes yielding F 3 = 0.23U 2 , and 
0.33(7 3/2 , respectively. 



0.2 

o.i 







-0.1 



-0.2 

0.2 0.4 0.6 jj 

FIG. 6. Parametric dependence of the three-particle spec- 
trum as a function of U, around the band center. The change 
in spectral correlation can be seen. For U < 0.2, many levels 
seem to cross (weak level repulsion), while for larger U, the 
crossings are avoided (strong level repulsion). 




In the numerical data, one observes indeed a transition 
at U ~ 2 to a regime in which cx U 3 ^ 2 with a pref- 
actor whose order of magnitude coincides again with the 
expected value. 

The upper points in Fig. |^ show the participation ra- 
tio R = {YZ=\ |C/3«| 4 ) -1 . The behavior of R is quite 
similar to the one of the spread width T3. The par- 
ticipation ratio, which goes to R = 1 at U = 0, in- 
creases proportional to the square of the interaction in 
the regime 0.25 < U < 1.5 and, as T3, exhibits the 
signature of the ergodic regime for U > 1.5. When 
calculating the participation ratio R, one has also to 
take into account, that not all of the states which are 
in the available energy interval can participate. Again, 
one has to replace the level spacing A3 by to ob- 
tain R « 7iT 3 /Af m 2Ti 2 U 2 Z(M 3 (Af) 2 )- 1 « 28U 2 . 
In the ergodic regime, the self-consistent estimate pre- 
sented above gives R « 212J7 3 / 2 . However, at low U, a 
difference arises since R = 1 and = at U = 0. 

The ratio -R/T3 w 260 in the quadratic regime is much 
smaller than the one expected from a democratic partic- 
ipation of the Fock states according to (||), which yielded 
R/T 3 = 7r/A 3 Ri 850. This is a consequence of the 
fact that individual states can show strong fluctuations 
around ([5]), thereby lowering the participation ratio. Fur- 
thermore, the sparse structure of the Hamiltonian, and 
the resulting hierarchical structure of the eigenfunctions 
reduces the number of participating basis states as has 
been seen in Fig. 0. 



C. Spectral statistics 

The evolution of some energy levels in the center of the 
band as a function of U is shown in Fig. ^. First of all, 
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FIG. 7. Top: Level spacing distribution P(s) for three 
spinless Fermions at U = 0.1, U = 0.3 and at U = 1. For 
comparison, the Poisson and GOE distributions are plotted 
also. Bottom: The integrated level spacing distribution / as 
a function of the interaction for two and three particles. 

the positive slope of all of the levels is visible. This is due 
to the diagonal matrix elements of the interaction which 
lead to a shift in energy of the order of 0.45L 7 as expected 
from the typical size of these elements. 

But we can also observe changes in the statistical be- 
havior of the spectrum. At low U, there are strong fluc- 
tuations in the level spacing while the spectrum becomes 
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FIG. 8. T,2(E) for the spectrum of two (top) and three 
(bottom) spinless Fermions in a 3 x 3 x 3 disordered lattice 
(V = 1, W = 2) for different values of U. 

more rigid around U « 0.4. This happens when neigh- 
boring states in energy (which usually are not directly 
coupled by the interaction) become correlated. 

This Poisson-Wigner transition can be studied more 
systematically by calculating the local fluctuation prop- 
erties of the spectrum. The first quantity we look at is 
the level spacing distribution P(s) in the center of the 
many-body spectrum around E = 0. P{s) for three par- 
ticles is shown for a few values of U in Fig. [?] (left). A is 
the corresponding mean level spacing (A2 or A3). At low 
U, P(s) is quite close to the Poisson limit of uncorrelated 
levels, while it tends towards the universal Wigner result 
(GOE) at stronger U. This transition is described in a 
quantitative way by the integral 

/:= / dsP(s), (24) 
J 1/2 

which is shown in the right hand side picture of Fig. |7]. 
The transition is much more abrupt for TV = 3 than for 
N = 2. This agrees with the expected (7-dependence of 
the Rabi frequency (\U\ for TV = 2 and \U\ 2 for TV = 3) 
which characterizes the range of the level repulsion for 
weak U . 

The energy scale Ejj characterizing the universal spec- 
tral Wigner-Dyson rigidity can be extracted from the 
variance T, 2 (E) =< N(E) 2 > - < N(E) > 2 of the num- 
ber of energy levels in an interval of width E. Comparing 
the behavior of the spectrum-to the GOE-behavior, one 
can identify this energy EuV3 (up to which the GOE- 
rigidity can be observed and above which one can see 




FIG. 9. The characteristic energy E\j up to which the 
2/3-particle spectrum exhibits the universal Wigner-Dyson 
rigidity. 

significant deviations). ^(E) is shown in Fig. ||, near 
the band center. One can observe, for two as well as 
for three particles, the Wigner-Dyson rigidity up to the 
energy scale Ejj where the spectrum becomes less rigid. 

For TV > 3 particles, one can see that the characteristic 
energy Ejj does not coincide with the spread width T 3 
of the eigenfuacjtions (as it is the case for two particles 
when T 2 > A 2 t3). 

For TV — 3, the crossover from the Poissonian behavior 
of uncorrelated levels to the universal behavior of the 
GOE is sharper. As can be seen in Fig. ^, Ey increases 
as U 2 for weak interaction. This corresponds to Rabi- 
oscillations due to second order coupling between nearby 
Fock states. This |?7 2 |-increase when Ejj < A3 for TV = 3 
is the analog of the \U\ -increase observed for TV = 2 when 
Ejj < A2. Above U c i, Ejj seems to linearly increase as a 
function of U. 



D. Large {/-limit 

For large U and small M, it is instructive to consider 
the Fock basis built of the on-site orbitals. In this basis 
the 3-particle states, and the role of U, can be classified 
according to the number of next-neighbor configurations. 
For a 3 x 3 x 3 cubic lattice, we have 1746 out of the 2925 
basis states where no next neighbor pairs occur and the 
Fock state is not shifted when the interaction increases. 
For 1008 basis states, there is one pair of particles nearby 
and the energy of the Fock states is E(U — 0) + U. For 
the remaining 171 Fock states the particles are clustered 
such that their energy increases like E(U = 0)+2U. Since 
the non-diagonal elements in this representation are only 
due to one-particle kinetic energy which does not depend 
on U, the different shifts of the Fock states lead to a 
splitting of the band into three parts, 1746 states around 
E = 0, 1008 states around E — U and 171 states around 
E = 2U. 

This can be seen for U — 15 in Fig.|l(], where the 
integrated density of states I DOS = J_ dE'p(E') is 
plotted for one realization of the disorder and different 
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result for Ejj w T\ in this non-ergodic Wigner-Dyson 
regime, one gets for three particles U c i — U c i ~ 0.8 since 
Af w 3.3A 3 . 

In summary, we find a clear evidence of the existence 
of two distinct thresholds U c i and U C 2, as well as a strong 
indication that the self-energy of an individual Fock state 
cannot be evaluated by perturbation theory for gjy > 1. 



L — ^ 1 1 1 1 — 1 

-10 10 20 30 „ 

E 

FIG. 10. The integrated density of the three-particle 
states for different values of the interaction. 

values of the interactions. The density of states can be 
fitted by a sum of three Gaussians centered at E = 0, 
U, and 2U with weights corresponding to the above 
mentioned numbers of occurrence of these shifts when 
U w 15. For U ~ 1, the interaction induced gaps in the 
spectrum are totally removed by the one-particle kinetic 
contributions. 



E. Evidence for the existence of two thresholds 

The numerical results for the three-particle spectrum 
clearly exhibit two different characteristic interaction 
strengths. First, the local fluctuations of the level spac- 
ing changes from Poisson to GOE when U ~ 0.3, as it can 
be seen from Fig. 0. This is consistent with the results 
obtained for E2 where Ejj « A3 when U » 0.3. Thus, 
adjacent energy levels are correlated when U > 0.3. This 
gives lid ~ 0.3. Below U c i, Ejj increases as U 2 in agree- 
ment with Rabi-oscillations in the perturbative regime. 
Above U = U c i, 33 = E1//A3 > 1 and the perturba- 
tion theory breaks down. An indication for this is that 
Ejj oc U and not Ey oc cx U A , as implied by the 
perturbation theory. Though we have no explanation for 
this linear behavior, we emphasize that this absence of a 
{/ 4 -behavior gives a strong hint that the range of validity 
of perturbation theory is limited to U typ < A| ff . 

Note that nothing striking is observed in the behavior 
of the total spread width T 3 of the states at U = U c \- 
There seems to be a crossover from R oc U to R oc U 2 
in the participation ratio R, but since R = 1 at U = 0, 
it might also be the signature of a saturation which is 
not connected to a characteristic energy scale but which 
becomes irrelevant when R 3> 1. 

However, there is something interesting happening at 
U ~ 1.5 where both, r 3 and R undergo a transition 
from a regime where they increase as U 2 to a regime 
where they increase as U 3 ^ 2 . This is clear evidence for 
a change in the structure of the eigenstates as it is ex- 
pected from our theoretical considerations at U = U C 2- 
At this strength of the interaction, we expect the indirect 
spread width to be of the order of ASp . The problem 
is now to estimate T\ when U > U c \ since the pertur- 
bation theory does not work. If one uses the numerical 



IX. IMPLICATION FOR QUANTUM 
LOCALIZATION IN REAL SPACE 

Very recently, scaling-type concepts have been ap- 
plied to two .particles with a local interaction. When 
Shepelyanskyt-3 had pointed out that in insulating sys- 
tems certain two-particle wave-functions could be delo- 
calized with respect to the one-particle states, Imry ex- 
tended the Thouless Block-scaling pictured and intro- 
duced a "two-particle conductance" 32 = Tu/A2, where 
T\j is the decay rate of the states in boxes of size L\ 
due to the interactive coupling to other boxes. This Ttj 
is identical to the T2 characterizing the spectral fluctu- 
ations for a given block. A2 is the mean spacing of the 
two-particle spectrum in a block. Using Fermi's golden 
rule for the estimate of Ttj, this yields a pair localization 
length L2 oc U 2 L 2 , in agreement with the results of Shc- 
pelyansky. The existence of this delocaii|a|tion effect has 
been confirmed in numerical studies" 6 H 17 1'B3. 

While in this paper we have considered the metallic 
phase, away from the localization transition, our treat- 
ment may form the basis for finding the localization 
length for an A-particle system (where N > 3). This 
localization length too is increased for large enough in- 
teractions. The scaling theory prescription for obtaining 
this localization length is straightforward, in principle. 
One has to increase the system size, L, and watch the 
L-dependence of both the effective spacing A^f of the 
A^-particle levels which are re-organized by the interac- 
tion when L > Li and the energy Ejj above which the 
Wigner-Dyson rigidity does not apply for this subset of 
levels. The L at which Ejj / A<ff becomes of order unity 
is the N-particle localization length. The parametric de- 
pendence embodied in our estimate of gs (<7at) can be 
used for this purpose. We expect the delocalization to 
become even stronger for N > 3. It is tempting to sug- 
gest that jsyph effects are at the origin of the recent ob- 
servationsE3c3 of a two dimensional metallic phase driven 
by the interactions in Si-MOSFET. 

X. QUASI-PARTICLE LIFETIME AND 
LOCALIZATION TRANSITION IN THE 
FOCK-SPACE 

We conclude by stressing the analogies and the dii. 
ferences between this study and the recently proposedEil 
approach to quasi-particlc lifetime in an isolated system. 
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In these two studies, two characteristic energies are iden- 
tified and a transition appears when their ratio is of or- 
der unity. This is the transition from weak perturbative 
mixing to the golden-rule decay, as mentioned in the in- 
troduction. This is also the threshold where the pertur- 
bation theory in U or e breaks down. In this sense, the 
two studies use very similar concepts, but the considered 
characteristic scales are not the same. This is because 
two different situations were considered (quasi-particles 
in a finite-density Fermi gas, vs. a dilute gas of "real" 
particles). It may be argued, as in RefEj, that the differ- 
ence between these two systems is mainly in the counting 
of the densities of excited states, and that in principle 
both could be treated by similar scaling considerations. 

We have mainly discussed what is the ratio (Thouless 
number) which controls the transition from Poisson to 
Wigner in the bulk of the TV-body spectrum. The thresh- 
old U Cl corresponds to Ut yp ~ A| a . Therefore, for very 
few particles where A>j ~ A2, and g n ~ gi ~ 1 at U c x, 
this first threshold is eventually related to the ratio g-i 
of the two-particle decay width T2 over the two particle 
spacing A2. This is not due to the fact that Ajv is not a 
relevant energy scale, but because this is the contribution 
of order P = N/2 or (N + l)/2 (depending on the parity 
of N) in the perturbation theory of the decay width Tjv 
which matters for the spectral fluctuations at the scale 
Ajy- Therefore, as previously proposed by Shepelyansky 
and Sushkov, we just need to have the two particle lev- 
els well coupled by the interaction in order to have the 
same thing for the iV-body levels. In this case the con- 
ditions for the establishment of Wigner-Dyson rigidity 
and to have an effective decay are similar and both are 

Utyp ~ A| . 

For the lifetime of quasi-particles in a zero-dimensional 
Fermi system, the relevant ratio is made from the two dif- 
ferent scales: the decay width of a single quasi-particle 
r sp (e) (resulting from the disintegration of a single quasi- 
particle into two quasi-electrons and a hole) and the ac- 
cessible three quasi-particle level spacing Aa(e). In this 
connection, a quasi-particle is considered at an excita- 
tion energy e above the Fermi energy of an isolated sys- 
tem. The corresponding Fermi vacuum is assumed to 
always provide a new electron-hole pair at each inter- 
action process, such that the relevant decay width is 
not the one characterizing the decay of a certain Slater 
determinant to those with the same quasi-particle con- 
tent, but to those with a quasi-particle content increased 
by a new electron-hole pair. It should be emphasized 
that this does not correspond to the dilute limit dis- 
cussed in our work, but to a limit of finite density of 
particles, such that the Fermi vacuum can be consid- 
ered as an unlimited reservoir of particle-hole excitations. 
F sp (e) w C/ t 2 yp /A 3 (e) w A^Xf ) defines the (second) exci- 
tation threshold e* in Ref.U. The threshold is obtained 
when the typical magnitude of the interaction matrix el- 
ement U t yp is of order A 3 (e), unlike thejialue A^ s which 
was relevant in our dilute limit. In Ref.liil, it is suggested 



that the two problems of level statistics and golden rule 
decay are unrelated: that derealization in Fock space 
does not mean that the spectrum should have Wigner- 
Dyson statistics. For a Cayley tree, this disagrees with 
the results of recent supersymmetric calculationsc9 using 
a non- linear sigma model formulation. For the "dilute" 
case of a small number of particles in a large volume, it 
follows from our work that these two properties are very 
intimately related. 

In this work, the transition^] from an uncorrelated 
spectrum to a fully ergodic one seems to occur in two 
stages: U c \ and U C 2- We do not have a qualitative un- 
derstanding of the behavior in the intermediate regime. 
Whether this regime will survive with increasing number 
of particles, and what is the precise r elationshipp] with 
Ref.EJ are questions which deserve further investigation. 
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